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Abstract 

We study the electronic structure of graphene in the presence of 
either sevenfolds or eightfolds by using a gauge field-theory model. 
The graphene sheet with topological defects is considered as a nega- 
tive cone surface with infinite Gaussian curvature at the center. The 
density of electronic states is calculated for a single seven- and eight- 
fold as well as for a pair of sevenfolds with different morphology. The 
density of states at the Fermi energy is found to be zero in all cases 
except two sevenfolds with translational factor M ^ 0. 
PACS: 73.22-f,73.22-Dj 

Nowadays graphene and related graphene-based nanostructures are the sub- 
ject of great interest. It is well known that some number of fivefolds must be 
included in the honeycomb graphene lattice to produce various nanostruc- 
tures with positive Gaussian curvature of the surface: nanocones, nanohorns, 
closed nanotubes, fullerenes, etc. There is another type of possible defects, 
sevenfolds and eightfolds, which results in a surface with negative Gaussian 
curvature. Experimentally the presence of sevenfolds in various nanostruc- 
tures was observed by Iijima []]. The structures with sevenfolds only are 
rather exotic [2], however, the presence of fivefold-sevenfold pairs is expected 
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in some cases and, what is important, leads to new features of the graphene- 
based materials. For example, in the carbon nanotube this pair changes chi- 
rality, which allows one to connect metallic and semiconducting nanotubes 
without breaking of 7r-bonds (see, e.g., [3]). Therefore the question of the 
electronic structure of carbon nanostructures containing sevenfolds has the 
definite theoretical and practical interest. 

Earlier computations on the carbon lattice jl] shows the general decrease 
of the density of states (DoS) near the sevenfolds in the sixfold plane. Notice 
that the influence of the curvature was not taken into account in [3]. Effects 
of both topological defects and local curvature on the electronic properties 
of planar graphene has been recently studied in [SI E] within the effective 
low-energy field-theory model proposed to study the electronic structure 
of slightly curved graphene sheets with an arbitrary number of pentagon- 
heptagon pairs. 

In the present paper, we investigate the electronic structure of the graphene 
plane with either a single sevenfold or two sevenf olds/one eightfold inserted 
(which corresponds to the negative 60° and 120° disclination, respectively). 
The surface of the plane is assumed to be free to bend but impossible to 
stretch. This results in appearance of the "surplus-angle cone" surface with 
infinite (negative) Gaussian curvature in the center, but with absence of the 
rotational symmetry. We use a field-theory model suggested in EJ [9] and, 
in final form, in [TO], which includes two gauge fields to take into account 
both the presence of sevenfolds and the influence of surface curvature. 

Let us consider the planar honeycomb lattice, which can be bent but can 
not be stretched. A single sevenfold can be inserted in this lattice by using 
the standard "cut-and-paste" procedure. The surface of the lattice due to 
the presence of negative disclination will appear to be curved (see Fig. 1). 
Let us assume the curved surface in the form 

r — (r cos ip sin 9((p),r sin sin 6 ((f), r cos 9 (<p)), (1) 

where 

%>) = f cos(2^) + |. 

After the cut-and-paste procedure, the length of the line of a constant dis- 
tance R from the sevenfold (which is also the line of intersection of the 
3d-sphere and the surface) is larger than for the same line in the planar case. 
Since the sevenfold leads to an additional 60° sector, the length should have 
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Figure 1: The curved surface of the honeycomb lattice with a single sevenfold 
inserted at the center (a "surplus-angle cone"). 



the value (7/6)2ttR. One can define it from the condition 

R ^ v/sin 2 9(tp) + (6'(<p)yd<p = 2nR 7 -, (2) 
Jo o 

which allows us to fix the only parameter a. The metrical tensor for the 
described surface is written as 

g„ = l,g w = r 2 [sin 2 6>(y?) + a 2 sin 2 (2y?)], g rip = 0, (3) 

and the non-zero metrical connection coefficients read 

r r w = -W.rSr = r ^ = Vr,r^ = (i/2)^in^ w . (4) 

The two-folds can be choosen in the form 

e r = 1> e % = a/5W>; e ]p = e r = 0- (^) 

By using the known relation (k> M ) = e l u D fl e l3l/ where an ordinary covari- 
ant derivative includes the metrical connection (J4]), one finds the spin 



3 



connection coefficients as 

K) 12 = -K) 21 = - 9 -f- (6) 

The model Dirac equation has the form (see [10J) 

- i<7 a e£(V„ - zaj - iWgV* = ^ fe , (7) 

where is measured from the Fermi energy, a a (a = 1,2) are the conventional 
Pauli matrices, a k (k = K, K_) and are two gauge fields (non-Abelian 
and Abelian, respectively), and = <9 M + (/x = r, <£>) with 

tyi = ^ 6 [(7 o ,0- 6 ]. 



In our case, one obtains 



2r 



To simplify the problem the gauge fields are taken in the same form as for 
the isotropic case: 



a 



±1/4, W v = -1/6, a k r = W r = 0. (8) 



The Dirac equation ([7]) can be diagonalized in respect to the K/K-- part 
(see [10] for detail), so that finally it takes the form 

- ia^dr + — )-0 ~ ia v ~ iWp)^ = ( 9 ) 

One can easily see that for a = (jSD reduces to the planar Dirac equation 
in the polar coordinates. Unlike the planar case, in g w is a function of 
the polar angle (p. This markedly complicates the analysis. Let us make the 
separation of variables in (Q assuming 

)<%>)■ 

The function $ should neglect the dependence of the wavefunction on (p in 
(Q. Therefore, we find the equation for $ in the form 

$' i 

— -ia v - iW v = i] Vsin 2 G(<p) + a 2 sin 2 (2<^), (10) 



where j is a constant. Introducing an auxiliary function 

ctp i 
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G(ip) = j ^sin 2 #(0) + a 2 sin 2 (20)d0, (11) 

one can write down $ in the form $ = exp^a^if + iW^ip) exp(ijG(ip)). Let 
us note that G(2it) = 27r(7/6) due to ([2]). The boundary condition for the 
spinor wavefunction has the form ip((p + 2tt) = —ip(ip). This allows us to find 
the constant j as 

j =3n = \ [1/2 + (n-a v - W v )\, n = 0, ±1, . . . (12) 

and 

<£>„ = exp(«(a ¥P + W-^V) 7= , (13) 

V 27T 

which finally leads to the system of coupled equations for u n ,v n in the fol- 
lowing form: 

~ i(dr + 7T + —) v n{r) = Eu n (r), 
2r r 

-Kdr + - -)u n {r) = Ev n (r). (14) 
2r r 

Notice that this system is equivalent to the planar Dirac equation with the 
only difference in the values of j n . The solution to (|T4l) is written as 



U, 



CJ u ^ 1/2l (\E\r), v n = ±iCJ ljn+1/2l (\E\r), (15) 



where J are the Bessel functions, and C is the normalization constant. The 
density of electronic states per unit area (LDoS) can be found from the 
properly normalized wavefunctions. The normalizing condition 

C 2 J r{u 2 + v 2 )dr = 1. 

allows us to find the normalization constant. Taking into account the asymp- 
totical behaviour of the solution (fT5l) ([^ + i/ 2 ( a; ) + Jn-i/2( x )] x — 2/7T when 
x — > oo) one can find the LDoS in the following form: 

LDoS(E,r) = — B - Y,(Mj n+ im(\m 2 + JeU n ~im(\E\r) 2 )- (16) 



Here Q (Q = 1/6 for a sevenfold) defines a surplus angle of a cone and 



j n = 6/7(1/2 + n + 1/6 ±1/4). 



(17) 



The summation runs over all values of n and e = ±1, which satisfies the 
condition ej n > 0. This restriction comes from the fact that the LDoS 
should be finite (or zero) in the limit E — > 0. In that way it overrides the 
normalization condition for the wavefunction at r — > 0. Notice that the LDoS 
does not depend on the angle coordinate tp despite the fact that the surface 
does not possess the axial symmetry. Accordingly, (fl6|) formally looks like 
the LDoS for a normal (deficit- angle) cone found in [TTJ. As a result, the 
only dimensional quantity in the LDoS is r while the energy can be reduced 
to the dimensionless parameter with respect to the energy unit hVp/r where 
Vp is the Fermi velocity. 

In the case of two closely situated sevenfolds or one eightfold the LDoS 
is also described by ffTBT) with = 1/3 and 



The factor M describes the morphology of defects: M = for an eightfold 
or when the shift vector (n, m) from one sevenfold to another satisfies the 
condition (n + m) mod 3 = 0, otherwise M = 1 (see [T0J [EE] for detail). The 
energy dependence of the LDoS near the center and far from the center is 
presented in Fig. 2. As is seen, the LDoS rapidly increases with energy near 
the sevenfold while far from the defect it has almost linear behavior. 

The LDoS as a function of the energy and the coordinate is shown in Fig. 
3. One can see that at the Fermi energy (E = 0) the LDoS is zero for both 
single defects while a finite value of LDoS is found for two sevenfolds with 
M ^ 0. This is provided by the term with j n = in the LDoS. The LDoS 
increases near the center (r pa 0) and becomes r-independent at far distances 
being linear in energy like in the case of planar graphene. 

To conclude, we have investigated the electronic structure of the graphene 
with one and two sevenfolds as well as with one eightfold in the framework 
of the field-theory approach. The curvature of the surface was taken into 
account by considering the "negative cone" with infinite and negative Gaus- 
sian curvature at the center. The calculations show that near defects the 
density of states is growing faster in comparison with the planar case. When 
the energy and/or a distance from the defect increases, the LDoS tends to be 
almost similar to the planar case. These findings are in general agreement 



j n = 6/8(1/2 + n + 1/3 ± (1/2 + M/3)). 



(18) 
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Figure 2: Fig.2. The local density of states per unit area (in the same scale), 
as a function of the energy near the center (top) and far from the center 
(bottom). 



with numerical calculations in [4J. At small distances from the defect an 
increase of the LDoS (per unit area) was found. Within our approach this 
increase can be explained by the influence of infinite curvature at the center 
(r = 0). In the case of two sevenfolds with the translational factor M^O 
the existence of the metallization near the defects was observed. For the 
eightfold, the metallization in the defect region is absent. This conclusion 
differs from the results of [I] where near the eightfold a non-zero value of the 
LDoS was found to vary (highly anisotropic) for different atoms on the same 
distance from the defect, thus forming local metal states. Such difference can 
be explained by two reasons: (a) the continual character of our model and 
(b) an isotropic approximation. Notice that in our study the LDoS does not 
depend on the angle <p due to excluded radial components of the fields jSJ). 
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On the other hand, the shape of the surface is propertly taken into account 
in our model, leading to curvature-induced peculiarities which are absent in 
the planar case [I]. 

This work has been supported by the Russian Foundation for Basic Re- 
search under grant No 08-02-01027. 
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Figure 3: The local density of states per unit area (in the same scale), as a 
function of the radial coordinate r and the energy E. The cases of a single 
sevenfold (top), single eightfold or two sevenfolds with M = (middle), and 
two sevenfolds with non-zero M (bottom) are shown. 



9 



